Quantum Zeno control of coherent dissociation 
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We study the effect of dephasing on the coherent dissociation dynamics of an atom-molecule Bose- 
Einstein condensate. We show that when phase-noise intensity is strong with respect to the inverse 
correlation time of the stimulated process, dissociation is suppressed via a Bose enhanced Quantum 
Zeno effect. This is complementary to the quantum zeno control of phase-diffusion in a bimodal 
condensate by symmetric noise (Phys. Rev. Lett. 100, 220403 (2008)) in that the controlled process 
here is p\ia&e- formation and the required decoherence mechanism for its suppression is purely phase 
noise. 

PACS numbers: 



Pair production processes are the key to the study 
of quantum correlations and to the preparation of non- 
classical states of light and matter. In optics, eminent 
correlated-photon experiments |H have led to the demon- 

\ stration of Einstein-Podolsky-Rozen (EPR) correlations 
In matter- wave optics, atom pairs may be produced 
by dissociation of ultracold diatomic molecules Q, po- 
tentially leading to interatomic quantum correlations [J] . 
In particular the dissociation of a Bose-Einstein Conden- 

' sate (BEC) of diatomic molecules into boson constituent 
atoms, is equivalent to optical parametric dovifnconver- 

\ sion in x^^-* nonlinear Kerr media Q, which is the lead- 
ing mechanism for the generation of correlated photons. 
The modulational instability implied within this analogy, 
may enable the generation of pair-correlated or number 
squeezed matter- wave beams j6j and open the way to the 
realization of Bose-stimulated superchemistry Q as well 
as to novel atom-interferometry techniques [8| . 

In this work we aim to study the interplay between 
coherent pair production and decoherence due to the 
coupling to an external bath. For matter-wave systems 
such dephasing mechanisms result from inelastic colli- 

\ sions between the Bose condensed atoms and particles 
(atoms or molecules) in the surrounding thermal cloud. 
It may also be introduced by random variation of the 
atom-molecule detuning via magnetic field fluctuations in 
Feshbach-resonance coupled systems [9] or via laser fre- 

\ quency fluctuations in optically coupled atom-molecule 

. BECs. 

Rather than study the detrimental effect of weak noise 
on the dissociation rate, we seek to deliberately use noise 
to demonstrate a Bose-stimulated quantum Zeno effect 
(QZE) [i3 resulting in the control of the atom-molecule 
system and the protection the coherent molecular BEC 
from stimulated dissociation. As the coherent dissocia- 
tion is sensitive to the relative phase between the atomic 
and the molecular modes, the continuous projection onto 
relative-number states (Fig. [Ija)) slows down the disso- 
ciation. The effect is complimentary to a recent study 
on the noise-suppression of phase-diffusion in double- well 
condensates However, whereas site- indiscriminate 

noise was used in (llj to protect against a process where 
phase-information is lost in favor of relative-number in- 
formation, here 'local' dephasing is utilized to arrest 
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FIG. 1: (Color online) (a) Noise-control of coherent dissocia- 
tion. Expectation values for SU(1,1) coherent states, are re- 
stricted to the upper- sheet hyperboloid {Kl) - {Kl) ~ {KD = 
. Phase-randomization projects onto the Kz number-states, 
thus slowing down the coherent dissociation of the initial 
atomic vacuum state, (b) Coherent (noise-free) atomic popu- 
lation dynamics, starting from the vacuum |0, 1/4) with x = 2 
(solid lines), x = 0.5 (dashed lines), and x = (dotted lines). 
Symbols correspond the linearized analytic expression of Eq. 
dll). Particle numbers are iV = 50 (o), 100 (□), and 200 (o). 

phase buildup by increased number fluctuations. 

We consider the parametric pair production Hamilto- 
nian 

H = ^a^a + ^{a)a)h + b'^aa), (1) 

where and w denote respectively, creation operators 
for atomic and molecular modes in the matter wave re- 
alization, g is a measure of the strength of the intermode 
coupling and e is the atomic binding energy in the ab- 
sence of coupling. The total atom number N = na + 2nb, 
where fia = a) a and fih = Wb, commutes with the Hamil- 
tonian ([T]), and is therefore conserved. 

The short-time dynamics of molecular dissociation 
retains SU(1,1) coherence [T2j. We define the three 
SU(1,1) generators, k+ = ataV2, K- = aa/2, ^ 
[K^,K+]/2 ~ + 1/4, with the complementary op- 

erators = {k+ + k^)/2,ky = {k+ - k^)/2i. The 
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mutual eigenstates of the Casimir operator and of 
Kz form the number basis set K'^\n, k) — k{k — l)|n, fc), 
k:,\n, k) = {k + n)\n, k), with n = 0, 1, . .■N/2, and the 
Bargmann index fc = 1/4 being the minimal value of 
when no atoms are present. The SU(1,1) coherent states 
\9, 6) are defined in relation to the number basis as 



1 - tanh"^ - 



k N/2 

E 



-tanh ( - ) exp(— i(/ 



(2) 



to suppress it should involve phase fluctuations which will 
project the state of the system into well-defined relative- 
number states. 

Such phase-noise may be implemented via stochastic 
variation of the atom-molecule detuning e or via colli- 
sions with thermal particles. Controlling the modula- 
tion rate or the density of thermal atoms thus provides a 
handle for controlling dissociation via the noise parame- 
ters. In the former case, we stochastically modulate the 
atom-molecule detuning (e.g. by introducing a fluctuat- 
ing magnetic field in a Feshbach resonance setup), thus 
modifying the system Hamiltonian as 



For such states the expectation values {Ki) = 
{9, (f>\Ki\9, (f>) are restricted to the upper-sheet of the hy- 
perboloid {K^f - (K^)^ - {Kyf = k^, shown in Fig^a). 

Expressing the Hamiltonian in the new operators, we 
get the simple form 
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(3) 



Starting from the atomic vacuum state \n = Q,k — 1/4) 
corresponding to an all-molecules BEC, it is therefore 
clear that the state of the system will remain an SU(1,1) 
coherent state as long as the molecular mode is macro- 
scopically populated and can be replaced by a classical 
c-number, so that the Hamiltonian ^ only generates 
Lorentzian boosts and rotations. 

Linearization about the atomic vacuum, assuming an 
undepleted molecular pump h — > N/2 (justified for 
short time dynamics where the depletion of the pump 
mode is negligible), results in the atom number evolu- 
tion 



Hs = H + C{t){K,--). 



(5) 



The stochastic term ^(i) is characterized by its zero mean 
(C(O)t = ^-iid the characteristic decay time {tc)noise 
of the noise correlation function {({t)C{t')). When 
itc)noise <C tc, wc Can assumc a delta-correlated white 
noise (c{t)C{t')) w r6{t - t'). The standard iteration 
of the Liouville equation then results in the Markovian 
kinetic master equation. 



(6) 



where p is the A'^-particle density matrix of the system, 
r is the noise intensity, and the Lindblad noise term C is 
the population difference. 
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t = 2y^h 



N -AK, + l. 



(7) 



where x = e/ {g\/2N) is the characteristic coupling pa- 
rameter and time has been rescaled as t = tg^N /2. 
Thus, as shown in Fig.[ljb), for weak atom- molecule cou- 
pling (x > 1) the coherent SU(1,1) dynamics correspond 
to small oscillations of the atomic population, whereas 
the near-resonance strong coupling (x < 1) evolution 
corresponds to the Bose-stimulated parametric amplifi- 
cation of emitted atom pairs Q . 

Regardless of the value of X: the short-time dynam- 
ics of Eq. Q corresponds to an initially quadratic de- 
cay with a vanishing rate as r — 0. The timescale on 
which the buildup of atomic population can be consid- 
ered quadratic, is the correlation time Tc = (1 — x^)~^^^ 
corresponding to tc = Tc/{gy^N/2) of order 1-10 ms for 
current BEC setups. This long correlation time enables 
the realization of a QZE [lO| which may slow down coher- 
ent molecular dissociation. Since the stimulated dissoci- 
ation process amounts to the formation of a well-defined 
phase between the atomic and molecular modes at the 
expense of increased number variance, the desired noise 



Taking the initial state to be |0, 1/4), the undepleted 
classical-pump approximation amounts to replacing b and 
with the fixed molecular population ^Jn/2. Using 
Eq. ([5]) in conjunction with the identities Tr(A[_B,C]) ~ 
Tr([A,B]C), Tr([^,S]) = 0, and the SU(1,1) commu- 
tation relations, we obtain dynamical equations for the 
expectation values Ki = Trjp-ftTi}: 



dr 



-2xKy - 16jK,, 
2xK, + 2K, - IQ-fKy, 

2Ky, 



(8) 



where 7 = F/ {g^N/2). Similarly, for the correlation 
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FIG. 2: (Color online) Atomic population dynamics starting 
from the atomic vacuum in the presence of noise. Panels 
(a),(b) are for subcritical x = 2 whereas panels (c),(d) are 
for supercritical x ~ 0.5. Dependence on noise-intensity is 
demonstrated in panels (a),(c) where N = 100 and 7 is 
(bold solid line), 5 (solid, o), 10 (dashed, A), and 20 (dash- 
dotted, o). Dependence on particle number is illustrated in 
panels (b),(d) where 7 = dynamics (bold lines) is compared 
with noise-suppressed dynamics with 7 = 10 (normal lines 
and symbols) for iV = 50 (solid, o), 100 (dashed. A), and 200 
(dash-dotted, o). Symbols correspond to the analytic QZE 
expression (|10|l . 



functions = TT{p{k^kj + KjKi)} - 2K,Kj we find. 



dr 

— A 

dT 



dr 
d 

d^ 
d 



-4xA,y - 327 (A,, - /^yy - 2K^) 



AX^.y + 4Ay, - 327 {^yy ^ A,, ^ 2 , 

2X (A^:, - Ayy) + 2A^2 - 647 {A^y + K^Ky) , 

-2xAy, + 2A:,j, - 167A^^, 

2xA,, + 2(Ayy + A,,)~167Ay, . (9) 



We note that Eqs. ([8|) and Eqs. dU are valid as long 
as the molecular field can be assumed to be undepleted 
{ua/N < 0.1) and do not make any further large-number 
assumption on the atomic operators Ki. They thus ap- 
ply also for the very early stages of dissociation when the 
atomic population comes purely from quantum sponta- 
neous emission. Unlike standard mean-field equations 
which have the atomic vacuum as an unstable stationary 
point, they accurately depict the amplification of and 
the parametric SU(1,1) phase-squeezing, when F = 0. 

If the frequent measurement condition F ^ t^^ is sat- 
isfied, Kx and Ky can be adiabatically eliminated in Eqs. 
([5]), in accordance with the QZE picture where coherence 
is not allowed to build up due to the introduction of noise. 
Substituting the initial values corresponding to the all- 
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FIG. 3: (Color online) The SU(1,1) Husimi distribution func- 
tion Q{0,(f)), projected onto the Kx,Ky plane. The initial 
state at r = (a) is a coherent vacuum, with equal ^J1TJ2 
variances in and Ky. In the absence of noise (7 = 0), 
a definite atom-molecule phase has been acquired at r = 3 
(b). By contrast, when noise is introduced with 7 = 0.1 (c) 
and 7 = 5 (d), squeezing is arrested and the initial coherent 
vacuum is protected. In all plots x = ^-nd A'' = 30. 



molecular state K = [0,0, 1/4], we find that the atomic 
population is then given by. 



fio_ _ ^ 

N ~ N 



exp 



7T 



XV4 + 167^ 



- 1 



(10) 



At the limit of on-resonance strong interaction, % — >■ 0, 
this expression reduces to the familiar QZE behavior: 



fia _ ^ 
N ^ N 



(11) 



corresponding to the slowing down of stimulated disso- 
ciation as the noise intensity F is increased. Interest- 
ingly, since the effective coupling strength g^/N is also 
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N dependent, the obtained QZE is Bose enhanced 
as discussed below. 

In Fig. [2] we show the early population dynamics ob- 
tained from direct numerical simulations of the dissocia- 
tion process according to Eq. (|ni). We show representative 
plots for X values above unity (panels a,b) and below it 
(panels c,d), corresponding to stable and unstable coher- 
ent dynamics, respectively. Compare first the noiseless 
evolution (solid bold lines) to the QZE-controUed evolu- 
tion in the presence of noise of varying intensity F (nor- 
mal lines), while keeping the particle number N fixed, 
as shown in Fig. [^h^ and Fig. Symbols depict the 
QZE prediction ([TU)) . showing excellent agreement with 
the exact numerical results. As expected, both stable os- 
cillations and unstable hyperbolic amplification are con- 
verted in the presence of noise, to similar F-controUed 
exponential buildup of atomic population, because the 
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FIG. 4: (Color online) Dynamics of the correlation functions 
Aij with A*' = 100 and x = starting from |0, 1/4), at various 
noise intensities: (a,b) 7 = 0, (c,d) 7 = 0.1, and (e,f) 7 = 5. 
Lines correspond to direct integration of the master equation 
((61), whereas symbols denote the undepleted-pump dynamics 
of Eqs. (O with initial conditions K = (0, 0, 1 /4) , A^x = 
Ayy — 1/4: , A^^ — A^y — A^z = Ay z = 0, corrBsp oudlug to 
the atomic vacuum. 



QZE dynamics is composed from repeated t < tc pieces 
which have a similar quadratic behavior in both cases. 

The dependence of the QZE on particle number is il- 
lustrated in Fig. [^b and Fig. where we plot the early 
dissociation dynamics for different values of N , fixing x 
and 7. The transition from a hyperbolic dissociation rate 
(|4]) which depends on N only through its onset time, to 
the QZE-suppressed dissociation which depends lin- 
early on N, introduces a factor of N/ \ogN between the 
coherent and noise-controlled processes. Thus the QZE 
is strongly Bose-aniplified. 

The effect of dephasing on the coherent dissociation 
dynamics, is shown in Fig. [31 using SU(1,1) Husimi quasi- 
probability distribution plots. The Husimi distribution 
is defined in terms of the SU(1,1) coherent states ([2]), as 
Q{e,(j)) = {e,(j)\p{t)\e,(j)) where p{t) is the state of the 
system at time t. The initial coherent vacuum state at 
t — (Fig. [Sja)) has equal variances in and Ky. 
Without noise, it evolves to an SU(1,1) squeezed coher- 
ent state (Fig. [31[b)). We note parenthetically, that un- 



like SU(2) or Glauber coherent states, SU(1,1) coherence 
and squeezing are not contradictory [l2j. Consequently 
a definite phase is formed between the atomic and molec- 
ular modes. In the presence of noise (Fig. [31[c),(d)) the 
squeezing is lowered significantly, with nearly complete 
protection of the initial molecular state available for suf- 
ficiently intense noise. This picture correlates well with 
the dynamics of fluctuations presented in Fig. 01 Without 
noise (Fig. |4l[a),(b)) we obtain the standard parametric 
amplification of the number fluctuations A^^, accompa- 
nied by an equal hyperbolic growth of Ky and Ayy at 
fixed Kx and A^x, and leading to decreased phase un- 
certainty (A(^)^ oc Axx/{Kx + Ky) as illustrated in 
Fig-EIb). In the presence of phase- noise the amplification 
of Ayy is suppressed and Axx begins to grow (Fig. jSJc) 
and Fig.|l^c),(d)) untfl at the QZE limit (Fig. jS^d) and 
Fig- Hie), (f)) both transverse fiuctuations Axx, ^yy grow 
symmetrically at a rate which decreases with F, suppress- 
ing the amplification of number fluctuations and the for- 
mation of a relative atom-molecule phase. We note that 
the undepleted-pump correlation equations @ agree well 
with the exact numerical simulations. 

To conclude, it is interesting to contrast our results 
here with previous work on the QZE suppression of 
phase-diffusion in a bimodal BEC_[ll| with underlying 
SU(2) spherical algebra. In Ref. [ll| we sought to slow 
down the dephasing between two BEC modes by tailoring 
noise which affects both site-modes symmetrically and 
projects the state of the system onto the initial coherent 
odd superposition of the two modes. By contrast, here we 
aim to slow down the process of phase-formation, using 
a much more tenable form of noise, which appears nat- 
urally via collisions or magnetic field fluctuations, and 
projects the state of the system onto the initial number 
state. As anticipated, we have observed that the coher- 
ent dissociation of a dimeric molecular BEC consisting of 
bosonic atoms can be prevented by such 'local' noise. As 
in the SU(2) case, the dissociation becomes less likely as 
the number of particles increases as N/logN. 

This work was supported by the Israel Science Founda- 
tion (Grants 582/07, 346/11) and by grant no. 2008141 
from the United States-Israel Binational Science Foun- 
dation (BSF). 



[1] A. Aspect, P. Grangier, and G. Roger, Phys. Rev. Lett. 
47, 460 (1981); A. Aspect, P. Grangier, and G. Roger, 
Phys. Rev. Lett. 49, 91 (1982); A. Aspect, J. Dalibard, 
and G. Roger, Phys. Rev. Lett. 49, 1804 (1982); D. C. 
Burnham and D. L. Weinberg, Phys. Rev. Lett. 25, 84 
(1970); C. K. Hong, Z. Y. Ou, and L. Mandel, Phys. 
Rev. Lett. 59, 2044 (1987); Z. Y. Ou and L. Mandel, 
Phys. Rev. Lett. 61, 50 (1988). 

[2] A. Einstein, B. Podolsky, and N. Rosen, Phys. Rev. 47, 
777 (1935). 

[3] T. Mukaiyama et al, Phys. Rev. Lett. 92, 180402 (2004); 
S. Diirr, T. Volz, and G. Rempe, Phys. Rev. A 70, 



031601 (R) (2004); M. Greiner, C. A. Regal, J. T. Stew- 
art, and D. S. Jin, Phys. Rev. Lett. 94, 110401 (2005); S. 
T. Thompson, E. Hodby, and C. E. Wieman, Phys. Rev. 
Lett. 94, 020401 (2005). 

[4] E. S. Fry, T. Walther, and S. Li, Phys. Rev. A 52, 
4381(1995); T. Opatrny and G. Kurizki, Phys. Rev. Lett. 
86, 3180 (2001); K. V. Kheruntsyan, M. K. Olsen, and 
P. D. Drummond, Phys. Rev. Lett. 95, 150405 (2005); B. 
Zhao et al, Phys. Rev. A 75, 042312 (2007); C. Gneiting 
and K. Hornberger, Phys. Rev. Lett. 101, 260503 (2008); 
N. Bar-Gill et al, Phys. Rev. Lett. 106, 120404 (2011). 

[5] P. D. Drummond, K. V. Kheruntsyan, and H. He, 



5 



Phys. Rev. Lett. 81 3055 (1998); J. Javanaincn and M. 
Mackic, Phys. Rev. A 59, R3186 (1999); A. Vardi, V. 
A. Yurovsky, and J. R. Anglin, Phys. rev. A 64, 063611 

(2001) ; M. W. Jack and H. Pu, Phys. rev. A 72, 063625 
(2005); U. V. Poulsen and K. Molmer, Phys. rev. A 63, 
023604 (2001); E. Rowen et al, Phys. Rev. A 72, 053633 
(2005); M. Ogren and K. V. Kheruntsyan, Phys. Rev. A 
82, 013641 (2010). 

[6] A. Vardi and M. G. Moore, Phys. Rev. Lett. 89, 090403 

(2002) ; K. V. Kheruntsyan and P. D. Drummond, Phys. 
Rev. A 66 031602(R) (2002); K. V. Kheruntsyan, Phys. 
Rev. A 71, 053609 (2005). 

[7] D. J. Heinzen, R. Wynar, P. D. Drummond, and K. V. 
Kheruntsyan, Phys. Rev. Lett 84, 5029 (2000); J. J. Hope 
and M. K. Olsen, Phys. Rev. Lett. 86, 3220 (2001); M. 
G. Moore and A. Vardi, Phys. Rev. Lett. 88, 160402 
(2002); I. Tikhonenkov and A. Vardi, Phys. Rev. Lett. 
98, 080403 (2007). 



[8] B. Yurkc, S. L. McCall, and J. R. Klaudcr, Phys. Rev. 

A 33, 4033 (1986); I. Tikhonenkov and A. Vardi, Phys. 

Rev. A 80, 051604(R) (2009). 
[9] R. A. Duine and H. T. C. Stoof, Phys. Rep. 396, 115 

(2004) ; E. Timmermans et al, Phys. Rep. 315, 199 
(1999); S. Brouard and J. Plata, Phys. Rev. A 72, 023620 

(2005) . 

[10] B. Misra and E. C. G. Sudarshan, J. Math. Phys. Sci. 

18, 756 (1977); L. A. Khalfin, JETP Lett. 8, 65 (1968); 

W. M. Itano, D. J. Heinzen, J. J. Bollinger, and D. J. 

Wineland, Phys. Rev. A. 41, 2295 (1990); A. G. Kofman 

and G. Kurizki, Nature (London) 405, 546 (2000). 
[11] Y. Khodorkovsky, G. Kurizki and A. Vardi, Phys. Rev. 

Lett. 100, 220403 (2008); Phys. Rev. A 80, 023609 

(2009). 

[12] L Tikhonenkov, E. Pazy, Y. B. Band and A. Vardi, Phys. 
Rev. A 77 063624 (2008). 



